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Standard Solar Models 

John N. BahcalP 

^■Institute for Advanced Study, Princeton, NJ 08540 

I review recent developments that affect standard solar model predictions of solar neutrino fluxes. 


1. INTRODUCTION 


A lot of progress has been made in understand¬ 
ing the robustness of solar model predictions since 
Neutrino 96 In this talk, I will first summa¬ 
rize the new ingredients and then give the current 
best-estimates. Then I will discuss the uncertain¬ 
ties in the predictions. 

Many of the results given here are adopted from 
the recent BP98 paper ^ . As we shall see in Sec¬ 
tion H, there is excellent agreement between stan¬ 
dard solar models calculated by different groups 
with different codes. 

If you want to obtain the numerical data that 
are discussed i n this talk, you can copy them from 
my Web site: http://www.sns.ias.edu/'^jnb . 


2. NEW INGREDIENTS 


In this section, I will first summarize the new 
and relevant results on nuclear fusion reactions 
and on the screening of nuclear reactions and 
then summarize the situation with respect to neu¬ 
trino cross sections. Finally, I will mention a few 
miscellaneous improvements that have been made 
since Neutrino 96. 


2.1. Nuclear reaction cross sections 

In January, 1997, the Institute for Nuclear 
Theory (INT) hosted a workshop devoted to de¬ 
termining the best estimates and the uncertain¬ 
ties in the most important solar fusion reactions. 
Thirty-nine experts in low energy nuclear ex¬ 
periments and theory, representing many differ¬ 
ent research groups and points of view, partic¬ 
ipated in the workshop and evaluated the ex¬ 
isting experimental data and theoretical calcula¬ 
tions. Their conclusions have been summarized in 
a detailed article authored jointly by the partic¬ 


ipants and published by the Reviews of Modern 
Physics [^. In general outline, the conclusions of 
the INT workshop paper confirmed and strength¬ 
ened previous standard analyses of nuclear fusion 
rates, although in a few important cases (for the 
^He(a, yj'^Be, ^Be(p, 7 )®B, and ^'^N(p, re¬ 

actions) the estimated uncertainties were deter¬ 
mined to be larger than previously believed. 

The largest change from what was used in 
the results presented at Neutrino 96 is the lower 
’’Be(p, 7 )®B cross section adopted by Adelberger 
et al. Previously, most authors construct¬ 
ing standard solar models used the Caltech (CIT) 
value for the ®B production cross section ||]. The 
difference between the INT and the CIT estimates 
of the ®B production cross section is due almost 
entirely to the decision by the INT group to base 
their estimate on only one (the best documented) 
of the six experiments analyzed by the CIT col¬ 
laboration. 

As we go alone, I will indicate how the princi¬ 
pal predictions of solar models depend upon the 
assumed ®B production cross section. 

2.2. Screening of nuclear reactions 

In one respect, the calculation of neutrino 
fluxes has simplified from Neutrino 96 to Neu¬ 
trino 98. The rather complicated expressions in 
the literature for the screening of nuclear fusion 
reactions by electrons and ions have been replaced 
by a simple analytic expression that was origi¬ 
nally derived by Salpeter for the case of “weak 
screening” only. 

Gruzinov and Bahcall employed a mean field 
formalism to calculate the electron density of the 
screening cloud using the appropriate density ma¬ 
trix equation of quantum statistical mechanics. 
Because of well understood physical effects that 
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are included for the first time in this treatment, 
the calculated enhancement of reaction rates does 
not agree with the frequently used interpolation 
formulae. For the sun, screening effects cause 
only small uncertainties in the predicted neutrino 
fluxes if the appropriate Salpeter formula is used. 

2.3. Neutrino cross sections and energy 
spectra 

Improved cross sections for neutrino absorption 
by gallium and by chlorine are now available as 
well as somewhat more precise standard (undis¬ 
torted) neutrino energy spectra. I summarize the 
results here. 

A number of authors do not use the best avail¬ 
able data for the neutrino cross sections and en¬ 
ergy spectra and some authors even give event 
rates (or SNU values) but do not say which cross 
sections and energy spectra they use so that one 
cannot interpret their results precisely. 

I have calculated [Q neutrino absorption cross 
sections for ^^Ga for all solar neutrino sources 
with standard energy spectra, and for laboratory 
sources of ®^Cr and ^^Ar; the calculations in¬ 
cluding, where appropriate, the thermal energy 
of fusing solar ions and use improved nuclear 
and atomic data. The ratio, R, of measured (in 
GALLEX and SAGE) to calculated ^^Cr capture 
rate is i? = 0.95 ± 0.07 (exp) -I- ^Qgg (theory) 
and was discussed extensively at Neutrino 98 by 
Gavrin and by Kirsten. 

I also calculated cross sections for specific neu¬ 
trino energies chosen so that a spline fit deter¬ 
mines accurately the event rates in a gallium de¬ 
tector even if new physics changes the energy 
spectrum of solar neutrinos. In order to make 
possible more precise analyses of event rates for 
neutrino scenarios which change the shape of 
the neutrino energy spectra from individual neu¬ 
trino sources, I evaluated and presented, for the 
first time, theoretical uncertainties for absorption 
cross sections at specific energies, as well as for 
the standard (undistorted) neutrino energy spec¬ 
tra. Also for use by people doing neutrino oscil¬ 
lation calculations, I calculated standard energy 
spectra for pp and CNO neutrino sources and pre¬ 
sented the results in Appendices and on my Web 
site. 


I note in passing that neutrino fluxes predicted 
by standard solar models, corrected for diffusion, 
have been in the range 120 SNU to 141 SNU since 
1968 §. 

A group of us have recently redetermined the 
standard shape for the ®B neutrino energy spec¬ 
trum 1^. The available data all seem to be con¬ 
sistent with each other within rather small un¬ 
certainties, so we were able to determine not only 
a best-fit energy spectrum but also two extreme 
spectra that are different by what we estimate is 
effectively ±3 ct. The uncertainties include esti¬ 
mates of the radiative and forbidden corrections. 
This improved spectrum yields a slightly different 
®B absorption cross section Q. 

I am somewhat nervous about the standard 
(undistorted) ®Be neutrino energy spectrum since 
it does depend upon rather old data Q . It would 
be very good if the a-particle energy spectrum 
from ®Be decay could be remeasured accurately 
in a new laboratory experiment. 

The situation for neutrino-electron scattering is 
good; cross sections are available that include 
electroweak radiative corrections. 

2.4. Miscellaneous improvements 

The standard solar model, BP98, that is dis¬ 
cussed in this report includes somewhat im¬ 
proved radiative opacities calculated by the Liv¬ 
ermore National Laboratory group, the so-called 
OPAL96 opacitie s [To[ |, and the improved OPAL 
equation of state |11[ . 

One improvement that I am rather proud of is 
new publicly available software that I have made 
available on my Web site (under Neutrino Soft¬ 
ware and Data) is a program to calculate so¬ 
lar neutrino rates and uncertainties. The prob¬ 
lem of calculating the uncertainties in the pre¬ 
dicted fluxes is somewhat complicated, especially 
since the uncertainties are asymmetric for some 
of the important input parameters. I decided to 
make this software publicly available (and there¬ 
fore polished it significantly) so that people could 
see explicitly what uncertainties were included for 
each parameter and how the uncertainties were 
combined. 

I have also polished my nuclear energy gener¬ 
ation code. The changes made in this code, al- 
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Table 1 

Standard Model Predictions (BP98): solar neu¬ 
trino fluxes and neutrino capture rates, with Icr 
uncertainties from all sources (combined quadrat- 
ically). 
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though they probably took me altogether a few 
weeks of programming and debugging time, ended 
up not changing neutrino flux predictions by more 
than a percent. 

3. BEST-ESTIMATE FLUXES 
AND EVENT RATES 

Table 1 gives the neutrino fluxes and their 
uncertainties for our best standard solar model, 
hereafter BP98. As discussed in the previous 
section, the solar model makes use of the INT 
nuclear reaction rates, recent (1996) Livermore 
OPAL radiative opacities, the OPAL equation of 
state, and electron and ion screening as deter¬ 
mined by the recent density matrix calculation. 

Figure ^ displays the calculated ’’Be and ®B 
neutrino fluxes for all 19 standard solar models 
with which we are familiar which have been pub¬ 
lished in the last 10 years in refereed science jour¬ 
nals. The fluxes are normalized by dividing each 
published value by the flux from the BP98 solar 
model H]; the abscissa is the normalized ®B flux 
and the ordinate is the normalized ’Be neutrino 
flux. The rectangular box shows the estimated 
3 ct uncertainties in the predictions of the BP98 


solar model. The abbreviations, which indicate 
references to individual models, are identified in 
the caption of Figure |^. 

All of the solar model results from different 
groups fall within the estimated 3cr uncertain¬ 
ties in the BP98 analysis (with the exception 
of the Dar-Shaviv model whose results have not 
been reproduced by other groups). This agree¬ 
ment demonstrates the robustness of the predic¬ 
tions since the calculations use different computer 
codes (which achieve varying degrees of precision) 
and involve a variety of choices for the nuclear 
parameters, the equation of state, the stellar ra¬ 
diative opacity, the initial heavy element abun¬ 
dances, and the physical processes that are in¬ 
cluded. 

The largest contributions to the dispersion in 
values in Figure are due to the choice of the nor¬ 
malization for S'i 7 (the production cross-section 
factor for ®B neutrinos) and the inclusion, or non¬ 
inclusion, of element diffusion in the stellar evo¬ 
lution codes. The effect in the plane of Fig. ^ of 
the normalization of Sn is shown by the differ¬ 
ence between the point for BP98 (1.0,1.0), which 
was computed using the INT normalization, and 
the point at (1.18,1.0) which corresponds to the 
BP98 result with the CIT normalization. 

Helioseismological observations have shown ||, 
0 that diffusion is occurring and must be in¬ 
cluded in solar models, so that the most recent 
models shown in Fig. [I] now all include helium 
and heavy element diffusion. By comparing a 
large number of earlier models, it was shown that 
all published standard solar models give the same 
results for solar neutrino fluxes to an accuracy of 
better than 10% if the same input parameters and 
physical processes are included @ 0 - 

The theoretical predictions in Table 1 disagree 
with the observed neutrino event rates, which 
are, see Ref. ||^ and the results presented at 
this conference by Lande, Gavrin, Kirsten, and 
Suzuki: 2.56 ± 0.23 SNU (chlorine), 72.2 ± 5.6 
SNU (GALLEX and SAGE gallium experiments), 
and (2.44 ±0.10) x 10®cm“^s“’ (®B flux from Su- 
perKamiokande). 

Bahcall, Krastev, and Smirnov jl^ have com¬ 
pared the observed rates with the calculated, 
standard model values, combining quadratically 
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Figure 1. Predictions of standard solar models since 1988. The figure shows the predictions of 19 standard 
solar models in the plane defined by the ^Be and ®B neutrino fluxes. The abbreviations that are used 
in the figure to identify different solar models are defined in the bibliographical item, Ref. [^. We 
include all standard solar models with which we are familiar that were published in refereed journals in 
the decade 1988-1998. All of the fluxes are normalized to the predictions of the Bahcall-Pinsonneault 98 
solar model, BP98 |Q. The rectangular error box defines the 3cr error range of the BP98 fluxes. The 
best-fit ^Be neutrino flux is negative. At the 99% C.L., there is no solution with all positive neutrino 
fluxes if the fluxes of CNO neutrinos are arbitrarily set equal to zero. There is no solution at the 99.9% 
C.L. if the CNO neutrinos are fixed at their standard solar model values. All of the standard model 
solutions lie far from the best-fit solution, even far from the 3a contour. 


the theoretical solar model and experimental un¬ 
certainties, as well as the uncertainties in the 
neutrino cross sections. Since the GALLEX and 
SAGE experiments measure the same quantity, 
we treat the weighted average rate in gallium 
as one experimental number. We adopt the Su- 
perKamiokande measurement as the most precise 
direct determination of the higher-energy ®B neu¬ 
trino flux. 

Using the predicted fluxes from the BP98 
model, the for the fit to the three experimental 
rates (chlorine, gallium, and SuperKamiokande) 


is 

Xssm( 3 experimental rates) = 61 . (1) 

The result given in Eq. (|^), which is approxi¬ 
mately equivalent to a 20tT discrepancy, is a quan¬ 
titative expression of the fact that the standard 
model predictions do not fit the observed solar 
neutrino measurements. 

The principal differences between the results 
shown in Table 1 and the results presented in 
our last systematic publication of calculated so¬ 
lar neutrino fluxes [|^ is a 1.3a decrease in the 
®B neutrino flux and l.lcr decreases in the ^^Cl 
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and capture rates. These decreases are 

due mainly to the lower ^Be(p, cross section 
adopted by Adelberger et al. |^. If we use, as 
in our recent previous publications, the Caltech 
rciT) value for the ®B production cross section 
1^, then the ®B flux is (^(®B,CIT) = b.lftog x 
10® cm-2s-\ SNU, 

Cl, CIT 

and = 13lt? SNU, all of 

Gallium 

which are within ten percent of the Bahcall- 
Pinsonneault 1995 best-estimates. 

4. UNCERTAINTIES 

In this section, I will first discuss the formal un¬ 
certainties in the solar model flux calculations and 
then review the strong constraints that helioseis¬ 
mology places on perturbations of the standard 
solar model. 

4.1. Uncertainties in the flux calculations 

Table 2 summarizes the uncertainties in the 
most important solar neutrino fluxes and in the 
Cl and Ga event rates due to different nuclear fu¬ 
sion reactions (the first four entries), the heavy 
element to hydrogen mass ratio (Z/X), the ra¬ 
diative opacity, the solar luminosity, the assumed 
solar age, and the helium and heavy element dif¬ 
fusion coefficients. The -|-p reaction causes a 
0.2% uncertainty in the predicted pp flux and a 
0.1 SNU uncertainty in the Cl (Ga) event rates. 

The predicted event rates for the chlorine and 
gallium experiments use recent improved calcula¬ 
tions of neutrino absorption cross sections ^,||. 
The uncertainty in the prediction for the gallium 
rate is dominated by uncertainties in the neutrino 
absorption cross sections, -1-6.7 SNU (7% of the 
predicted rate) and —3.8 SNU (3% of the pre¬ 
dicted rate). The uncertainties in the chlorine ab¬ 
sorption cross sections cause an error, ±0.2 SNU 
(3% of the predicted rate), that is relatively small 
compared to other uncertainties in predicting the 
rate for this experiment. For non-standard neu¬ 
trino energy spectra that result from new neu¬ 
trino physics, the uncertainties in the predictions 
for currently favored solutions (which reduce the 
contributions from the least well-determined ®B 
neutrinos) will in general be less than the val¬ 


ues quoted here for standard spectra and must 
be calculated using the appropriate cross section 
uncertainty for each neutrino energy [^|j^. 

The nuclear fusion uncertainties in Table 2 
were taken from Adelberger et al. |3|, the neu¬ 
trino cross section uncertainties from the 

heavy element uncertainty was taken from helio- 
seismological measurements p^ , the luminosity 
and age uncertainties were adopted from BP95 
fl^ ], the Icr fractional uncertainty in the diffu¬ 
sion rate was taken to be 15% [^, which is sup¬ 
ported by helioseismological evidence 0 , and 
the opacity uncertainty was determined by com¬ 
paring the results of fluxes computed using the 
older Los Alamos opacities with fluxes computed 
using the modern Livermore opacities [ p^ . To 
include the effects of asymmetric errors, the now 
public-available code for calculating rates and un¬ 
certainties (see discussion in previous section) was 
run with different input uncertainties and the re¬ 
sults averaged. The software contains a descrip¬ 
tion of how each of the uncertainties listed in Ta¬ 
ble 2 were determined and used. 

The low energy cross section of the ^Be ± p re¬ 
action is the most important quantity that must 
be determined more accurately in order to de¬ 
crease the error in the predicted event rates in 
solar neutrino experiments. The ®B neutrino flux 
that is measured by the Kamiokande |Q , Super- 
Kamiokande 0, and SNO experiments is, 
in all standard solar model calculations, directly 
proportional to the ^Be±p cross section. If the Icr 
uncertainty in this cross section can be reduced 
by a factor of two to 5%, then it will no longer be 
the limiting uncertainty in predicting the crucial 
®B neutrino flux (cf. Table 2). 

4.2. How large an uncertainty does helio- 
seismology suggest? 

Could the solar model calculations be wrong by 
enough to explain the discrepancies between pre¬ 
dictions and measurements for solar neutrino ex¬ 
periments? Helioseismology, which confirms pre¬ 
dictions of the standard solar model to high pre¬ 
cision, suggests that the answer is probably “No.” 

Figure || shows the fractional differences be¬ 
tween the most accurate available sound speeds 
measured by helioseismology and sound 
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Table 2 

Average uncertainties in neutrino fluxes and event rates due to different input data. The flux uncertainties 
are expressed in fractions of the total flux and the event rate uncertainties are expressed in SNU. The 
^Be electron capture rate causes an uncertainty of ±2% jl^ that affects only the ^Be neutrino flux. The 
average fractional uncertainties for individual parameters are shown. 
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speeds calculated with our best solar model (with 
no free parameters). The horizontal line cor¬ 
responds to the hypothetical case in which the 
model predictions exactly match the observed 
values. The rms fractional difference between 
the calculated and the measured sound speeds 
is 1.1 X 10“^ for the entire region over which 
the sound speeds are measured, O.Obi?© < R < 
O.95i?0. In the solar core, 0.051?© < R < 0.251?© 
(in which about 95% of the solar energy and neu¬ 
trino flux is produced in a standard model), the 
rms fractional difference between measured and 
calculated sound speeds is 0.7 x 10“^. 

Helioseismological measurements also deter¬ 
mine two other parameters that help characterize 
the outer part of the sun (far from the inner re¬ 
gion in which neutrinos are produced): the depth 
of the solar convective zone (CZ), the region in 
the outer part of the sun that is fully convective, 
and the present-day surface abundance by mass 
of helium (l^urf)- The measured values, l?cz = 
(0.713 ± 0.001)l?o ||, and Tsurf = 0.249 ± 0.003 
[[T^ , are in satisfactory agreement with the val¬ 
ues predicted by the solar model BP98, namely, 
l?cz = 0.7141?©, and Kiurf = 0.243. However, we 
shall see below that precision measurements of 
the sound speed near the transition between the 
radiative interior (in which energy is transported 
by radiation) and the outer convective zone (in 


which energy is transported by convection) re¬ 
veal small discrepancies between the model pre¬ 
dictions and the observations in this region. 

If solar physics were responsible for the solar 
neutrino problems, how large would one expect 
the discrepancies to be between solar model pre¬ 
dictions and helioseismological observations? The 
characteristic size of the discrepancies can be es¬ 
timated using the results of the neutrino exper¬ 
iments and scaling laws for neutrino fluxes and 
sound speeds. 

All recently published solar models predict es¬ 
sentially the same fluxes from the fundamental 
pp and pep reactions (amounting to 72.4 SNU 
in gallium experiments, cf. Table 1), which are 
closely related to the solar luminosity. Comparing 
the measured gallium rates (reported at Neutrino 
98) and the standard predicted rate for the gal¬ 
lium experiments, the ^Be flux must be reduced 
by a factor N if the disagreement is not to ex¬ 
ceed n standard deviations, where N and n sat¬ 
isfy 72.4 -h (34.4)/7V = 72.2 -H na. For a la (3ct) 
disagreement, N = 6.1(2.05). Sound speeds scale 
like the square root of the local temperature di¬ 
vided by the mean molecular weight and the ^Be 
neutrino flux scales approximately as the 10th 
power of the temperature . Assuming that the 
temperature changes are dominant, agreement to 
within Icr would require fractional changes of or- 
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Figure 2. Predicted versus Measured Sound Speeds. This figure shows the excellent agreement between 
the calculated (solar model BP98, Model) and the measured (Sun) sound speeds, a fractional difference 
of 0.001 rms for all speeds measured between O.O5i?0 and 0.957?©. The vertical scale is chosen so as to 
emphasize that the fractional error is much smaller than generic changes in the model, 0.03 to 0.08, that 
might significantly affect the solar neutrino predictions. 
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der 0.09 in sound speeds {3a could be reached 
with 0.04 changes), if all model changes were in 
the temperature^. This argument is conservative 
because it ignores the contributions from the ®B 
and CNO neutrinos which contribute to the ob¬ 
served counting rate (cf. Table 1) and which, if 
included, would require an even larger reduction 
of the ^Be flux. 

I have chosen the vertical scale in Fig. 1 to 
be appropriate for fractional differences between 
measured and predicted sound speeds that are of 

have used in this calculation the GALLEX and SAGE 
measured rates reported by Kirsten and Gavrin at Neu¬ 
trino 98. The experimental rates used in BP98 were not 
as precise and therefore resulted in slightly less stringent 
constraints than those imposed here. In BP98, we found 
that agreement to within Icr with the then available ex¬ 
perimental numbers would require fractional changes of 
order 0.08 in sound speeds (3 (j could be reached with 0.03 
changes.) 


order 0.04 to 0.09 and that might therefore af¬ 
fect solar neutrino calculations. Fig. 1 shows that 
the characteristic agreement between solar model 
predictions and helioseismological measurements 
is more than a factor of 30 better than would be 
expected if there were a solar model explanation 
of the solar neutrino problems. 

5. DISCUSSION AND CONCLUSION 

Three decades of refining the input data and 
the solar model calculations has led to a predicted 
standard model event rate for the chlorine exper¬ 
iment, 7.7 SNU, which is very close to the best- 
estimate value obtained in 1968 |^, which was 
7.5 SNU. The situation regarding solar neutrinos 
is, however, completely different now, thirty years 
later. Four experiments have confirmed the detec¬ 
tion of solar neutrinos. Helioseismological mea- 








surements show (cf. Fig. 1) that hypothetical de¬ 
viations from the standard solar model that seem 
to be required by simple scaling laws to ht just 
the gallium solar neutrino results are at least a 
factor of 40 larger than the rms disagreement be¬ 
tween the standard solar model predictions and 
the helioseismological observations. This conclu¬ 
sion does not make use of the strong evidence 
which points in the same direction from the chlo¬ 
rine, Kamiokande, and SuperKamiokande exper¬ 
iments. 

The improvement in helioseismological mea¬ 
surements over the past two years, from Neutrino 
96 to Neutrino 98 (cf. Figure 2 of the Neutrino 96 
talk Q with Figure of this talk), has resulted 
in a five-fold improvement in the agreement with 
the calculated standard solar model sound speeds 
and the measured solar velocities! I believe that 
this improved agreement is yet another reason to 
believe that standard solar models reliably pre¬ 
dict solar neutrino fluxes. 

I am grateful to Y. Suzuki for special efforts 
that made possible my attendance at Neutrino 
98 and which enriched the scientific experience 
of this extraordinarily exciting conference. I ac¬ 
knowledge support from NSF grant ^PHY95- 
13835. 
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